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We predict an unconventional magnetic ground state in AyFei,6^e2 with x Fe-vacancy 
superstructure under hydrauUc external pressure based on first-principles simulations. While the 
Fe-vacancy ordering persists up to at least ^ 12GPa, the magnetic ground state goes at ^lOGPa 
from the BS-AFM phase to a Neel-FM phase, a ferromagnetic arrangement of a " Neel cluster^\ 
The new magnetic phase is metallic, and the BS-AFM to Neel-FM phase transition is accompanied 
by a sizable structural change. The two distinct magnetic phases can be understood within the 
extended J1-J2 Heisenberg model by assuming a pressure-tuned competition between the intrablock 
and interblock nearest-neighbor couplings of iron moments. 

PACS numbers: 75.25.-j,71.20.-b,75.10.Hk 



Introduction. The anti-PbO-type FeSe is the sim- 
plest Fe-based superconductor with Tc ^ 8K[i|. Re- 
cently, a class of new iron superconductors AyFe2-x^^2 
(A=Tl,K,Rb,Cs) with enhanced Tc - SOKjl-g has at- 
tracted intensive interest. These materials are struc- 
turally similar to the 122-type iron-pnictides[5| with the 
FeSe layers intercalated by the A atoms, leaving cer- 
tain amount of Fe- vacancies in the Fe-square lattice. 
The iron-vacancies are expected to order in some peri- 
odic superstructures at certain x values [3]. Among the 
proposed superstructures, the a/5 x a/5 vacancy order- 
ing pattern(see in FIG. [1]), corresponding to x = 0.4, 
seems to be of special importance since it exists in most 
of AyFe2-x^^2 compounds as confirmed by experiments 
from neutron diffraction [6|-|9| to high-resolution transmis- 
sion electron microscope [10|. 

In addition to the a/5 x a/5 superstructure, a novel 
block-spin antiferromagnetism, the BS-AFM order shown 
in FIGdl^b) with large magnetic moment of irons, was 
discovered by the neutron diffraction experiments [6,]. A 
remarkable observation is the co-existence of supercon- 
ductivity and BS-AFM order in Ay¥e2-x^^2[^' It has 
been debated whether the co-existence is intrinsic at the 
microscopic level [11], or due to phase separation [12]. Pre- 
vious independent density functional theory (DFT) cal- 
culations show that the ground state has a BS-AFM order 
in the presence of the Fe-vacancy superstructure [l3|, Q . 
Since the magnetic ordering temperature, which is un- 
precedentedly high (Tn ~ 550i^), is close to the Fe- 
vacancy ordering temperature Ty ~ 580i^[6], it was also 
proposed that the Fe-vacancy ordering may be driven 
by magnetic exchange interactions so as to minimize the 
magnetic frustrations [l5|. 

All these interesting issues are closely related to the 
Fe-vacancy orderings and call for further experimental 
investigations. A key to resolve these issues is to clarify 
whether the BS-AFM is the only magnetic ground state 
or whether there are other magnetic ground states in the 



Fe-vacancy ordered compounds. Here we suggest to seek 
for these states by applying physical pressure, a clean pa- 
rameter to tune the lattice and electronic structures. We 
note a recent high pressure experiment on superconduct- 
ing sample of nominal K0.8Fe1.7Se2 compound, where the 
resistance hump is suppressed at a critical pressure of 8.7 
GPa and Tc is suppressed at a similar pressure [16]. The 
metallic phase in the high pressure regime was believed 
to arise from charge transfer between two different Fe va- 
cancy occupancies but its magnetic structure is unknown. 

In this paper, we study the lattice and magnetic struc- 
tures of AyYei^Q^e2 under high pressure up to 16 GPa 
by using first-principles simulations. We find that while 
the BS-AFM ground state (which is insulating at = 
0.8[l3,ll3) persists to lOGPa, a novel metallic magnetic 
phase, the Neel-FM phase, becomes the new ground state 
for higher pressures. Both the BS-AFM and Neel-FM 
phases can be described by the extended J1-J2 model 
while the magnetic phase transition is accompanied by a 
sizable structural change. 



(a) 




(b)BS-AFM 



FIG. 1: The geometry (a) and magnetic patterns for (b) BS- 
AFM (AFM2) and (c) Neel-FM (AFM4). The structural pa- 
rameters are as defined in FIGdl 



2 



Magnetic phases and ground state. For the first prin- 
ciple simulations, plane-wave basis and projected aug- 
mented wave methods are used as implemented in VASP 
code [I7|,ll8|. The Perdew, Burke, and Ernzerhoff flavour 
of generalized gradient approximation (GGA)[19] is em- 
ployed to calculate the electron exchange-correlation en- 
ergy. A 360 eV energy cutoff for the plane-waves and 
4x4x4 Monkhorst-Pack k-grid[20] are chosen to en- 
sure the total energy converges to 1 meV/Fe. The exter- 
nal pressure is introduced using the Pulay stress method. 
With the applied pressure, both the internal coordinates 
and the lattice constants are fully optimized for each cal- 
culated magnetic and non-magnetic configurations until 
the total force on each atom is < 0.01 eV/A. For the 
density of state (DOS) calculations, a dense 16 x 16 x 16 
F-centered k-grid, as well as the tetrahedra method are 
used to obtain accurate energy gap sizes. In the follow- 
ing text, we report and discuss detailed results for A=T\^ 
while the validity of all conclusions is checked for A=K 
and Rb. 

We first identify the true ground state by examin- 
ing the relative energies of several possible magnetic 
configurations [13|, |21i]. As shown in FIG. Efa), the rel- 
ative energetic order of different magnetic configurations 
vary drastically with the applied external pressure. The 
AFM2 pattern (the BS-AFM phase, FIGlT^b)) remains 
the ground state until around lOGPa, where the origi- 
nally second highest AFM4 (FIGHJc)) pattern takes over 
and eventually becomes the new ground state on the high 
pressure side. In the AFM4 phase, the four spins in each 
Fe4 block form a Neel order, while all block configura- 
tions are parallelly aligned, thus we denote it the Neel- 
FM phase for simplicity. The new ground state is at least 
15 meV/Fe lower than all other calculated phases at 12 
GPa, and is at least 20 meV/Fe lower at 16 GPa. 

Structure variation. A close examination reveals the 
details of the vacancy-induced structural distortion dur- 
ing the process, as shown in FIG. [2jc)-(e). From GPa 
to 10 GPa, the interblock nearest-neighbor (n.n.) Fe-Fe 
distance c^pe-Fe compressed from 2. 86 A to 2. 64 A, more 
severely than the intrablock n.n. Fe-Fe distance dpe-Fe 
(from 2.62 A to 2.50 A). Meanwhile, ah Fe-Se-Fe angles 
(a, a' ^ a" ^ /3, and P') are slightly reduced. The forma- 
tion of vacancy superstructures also causes the relative 
height of Fe atoms zpe (measured from the nearest Tl- 
layer) to be slightly different in two neighboring blocks. 
However, this distortion is negligible {z^^ — z^^ < 0.08A) 
in the BS-AFM state until 10 GPa. 

The magnetic phase transition then occurs, together 
with a structural change, identified by the first abrupt 
volume change around lOGPa as shown in FIGl2fb). 
Across the transition, dpe-Fe expands to 2.61 A, d'p^_p^ 
is further reduced to 2.57A, and /3 jumps from 94.35° 
to 104.61°. The lattice constant a expands from 8.17A 
to 8.32A, while c is severely compressed from 13.24A to 
12.41A. The difference in z-pe is enhanced from 0.07A 



to 0.26A. After the transition, from 10 GPa to 14 GPa, 
dpe-Fe and c^pe-Fe ^^^^ slightly drops, but all Fe-Se-Fe 
angles start to increase, and the difference in z^e keeps 
increasing as well. 

Another abrupt change in Fe-Fe distances (and Fe-Se- 
Fe angles) could be identified from 14 GPa to 16 GPa 
(FIGEfb)), as the lattice constant a abnormally expands 
from 8.27A to 8.40 A and c collapses from 11. 95 A to 
11. 03 A. A closer examination reveals that the difference 
in ZFe expands abruptly from 0.32 A at 14 GPa to 0.73 A 
at 16 GPa. This severe change implies that the back- 
bone of vacancy superstructure is becoming unstable and 
that the superstructure may break down under such high 
pressure. Thus we focus on the magnetic phase transition 
around 10 GPa in the following discussions. 

Extended J1-J2 Heisenberg model. In order to under- 
stand the physics behind the pressure-induced magnetic 
phase transition we fit the energetics of the magnetic 
configurations using the extended J1-J2 model [13] (FIG. 
[SJa)) where Ji and J2 (or J[ and J2 ) are the intra- 
block (or the interblock) n.n. and next-nearest-neighbor 
(n.n.n.) exchange couplings, respectively. J2 is always 
AFM, for pressure from GPa to 16 GPa; while Ji, J[ 
and J2, being FM at GPa, become AFM at approxi- 
mately 8 GPa, 6 GPa, and 11 GPa, respectively. All Js 
increase monotonically with increasing pressure, except 
for J2 which dwells around -20 meV until lOGPa. We 
notice that for ambient pressure the fitted values of Js 
are compatible with the low energy spin wave excitations 
probed in the recent experiment [221], while in the high 
pressure Neel-FM phase, the fitted values are compatible 
with the Monte Carlo result [23]. 

Variations of these exchange interactions, while reflect- 
ing complicated electronic structures that evolve with 
pressure, can be attributed to a combined effect of the 
electron correlation ([/), the Hund's coupling (J//), the 
crystal field splitting, as well as various short-ranged hop- 
ping integrals at the microscopic five-orbital Hubbard 
model level. As inferred from the previous DFT study 
for the iron pnictides[24], Jh plays an important role in 
mapping out the local magnetic interactions though U 
involved in GGA calculations may be not large. Some 
insights can be gained if we adopt the intuitive results 
obtained from perturbation or Hartree-Fock mean-field 
theory 0,11^. The n.n. and n.n.n. exchange couplings 
are then contributed from two virtual processes where 
Fe-3d electrons hop between two sites directly, or via p- 
orbitals of the Se- atoms [24]. The contribution from the 
direct exchange depends strongly on the inter-atomic dis- 
tances, while the contribution from the indirect exchange 
depends strongly on both the inter-atomic distance and 
the Fe-Se-Fe angle. 

On one hand, both contributions decrease monotoni- 
cally with Hund's coupling and become FM at relatively 
large Jh/U. On the other hand, for fixed Jh/U^ the 
second contribution could be ferromagnetic when the an- 
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FIG. 2: Variation of magnetic configuration energies and structural parameters with respect to external pressure, (a) Total 
energy per unit cell of different magnetic phases. To enhance the visibility, the total energies have been renormalized using 
E{P) = E{P) -kP + Eo with k = 2.800 eV/GPa and Eq = 120.163 eV. (b) Unit ceh volume of BS-AFM and Neel-FM phases, 
where the black solid line indicates the actual volume per unit cell, (c)-(e) Various structural parameters as defined in Fig. 



gle 7 = {tt — P)/2 formed by the Fe-Se bond and the Fe 
plane is larger than a certain threshold value varying in 
between 35° ~ 40° approximately, depending on the de- 
tails of materials |26|. From our structural analysis, the 
threshold values of 7 calculated by the structure parame- 
ters for the corresponding fitted Ji , J( , and J2 are about 
36° ~ 42°, which reasonably fall into the FM region given 
by the mean- field approximation[26|. Meanwhile, the 7 
angle associated with J2 is the smallest (about 35° or 
less) due to the structural distortion. It explains why 
Ji, J(, and J2 are initially FM and increase under pres- 
sure, while J2 is AFM even for ambient pressure. It is 
interesting to remark that across the phase transition, P 
angle abruptly changes from 94.35° to 104.61°, leading 
to a sign change in the second contribution, and thus 
J2 becomes AFM around 10 GPa. From 10 GPa to 14 
GPa, not only Fe-Fe distances are reduced, but the Fe- 
Se- Fe angles (or 7) are also increased (reduced), thus all 
Js rapidly increases. 

Once the two competing phases, namely the BS-AFM 
and Neel-FM phases, are identified as shown in FIGEfa), 
the magnetic phase transition point should be at Ji = 
J{/2, when the energy of the Neel-FM phase is equal to 
that of the BS-AFM phase. Using this criterion and the 
least-square-fitting of the Js, we determine the critical 
pressure to be 10.53GPa. Hence the phase transition 
around 10 GPa is mainly driven by competition between 
Ji and J(, while the role of J2 or J2 is to stabilize the 
BS-AFM phase or the Neel-FM phases in the respective 



low or high pressure regimes |27|. 

Band structure and density of states. The transition 
from the BS-AFM to Neel-FM phases drastically changes 
the electronic structure of the material. It was shown 
that 74o.8Fei.6Se2 is an AFM insulator with the band gap 
- 400-600 meV [il,|lj]. Here, the gap size Eg and the 
magnetic moment per Fe atom mpe of Tlo.8Fe1.6Se2 in 
the BS-AFM or the Neel-FM states are calculated and 
shown in FIG.^b). 

In the BS-AFM phase, when the pressure increases 
from to lOGPa, the intra- and inter-block Fe-Fe n.n. 
distances greatly reduce, which enhance the hopping be- 
tween the n.n. Fe atoms, resulting in a decrease of Eg 
from ~400 meV to ~100 meV. Meanwhile, mpe starts 
from 2.81 /i^ at OGPa and gradually reduces to 2.40 
/iB at lOGPa. The rapid suppression of Eg compared 
with the small variation of mpe reveals possible vacancy- 
enhanced Mott physics [28|-|30|. 

By contrary, we find that the Neel-FM phase is always 
metallic as shown in FIG. 3(b). We have also checked 
that the metallicity is robust until U ~ 4 eV by us- 
ing GGA-\-U calculations. Therefore, for y = 0.8, the 
magnetic phase transition is associated with an insulator- 
metal phase transition. In the Neel-FM phase, mpe sets 
off from 2.52 /i^ at OGPa, and is reduced to 1.66 /i^ at 
10 GPa. It reaches the lowest value of 1.27 /i^ at 14GPa, 
where the collapsed phase transition occurs and restores 
mpe to 1.40 jj^B at 16GPa. It is important to notice that 
although the Neel block ensures vanishing net moment 
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FIG. 3: The variation of electronic structure with respect 
to external pressure, (a) Exchange- coupling constants fitted 
using the extended J1-J2 model, and (b) band gap size Eg and 
magnetic moment per Fe mpe- In panel (b), the open/filled 
circles and squares denote mpe and Eg at the BS-AFM/Neel- 
FM phases, respectively. 

in the classic ordering configuration, there is no overall 
q- vector ensuring the time-reversal symmetry, since the 
blocks are ferromagnetically aligned. 

Discussions. The magnetic states considered in this 
paper are based on the classical Heisenberg model, and 
the mpe, which is about 2.5/iB-3.0/iB in the BS-AFM 
phase, should be considered as the static local magnetic 
moment. This is in agreement with the neutron diffrac- 
tion experiment where a large iron moment 3.31/iB/Fe 
(for A=K) was reported [6]. Therefore quantum fluctua- 
tions should be suppressed in the ordered phases due to 
the large magnetic moment. One expects that quantum 
fluctuations play a larger role when the magnetic tran- 
sition point is approached, but the two magnetic phases 
should be robust in the presence of the vacancy super- 
structure. The magnetic phase transition is likely of first 
order because it is accompanied with a sizable change in 
the lattice constants. An adequate account of quantum 
fluctuations needed for a thorough understanding of the 
transition calls for future studies. 

Our results shed a new light in understanding the 
available high pressure experiment [16] if the resistance 
hump observed in most of the superconducting samples 
of AyFe2-xS^2 is interpreted as due to a phase separation 
involving a Fe- vacancy disordered superconducting phase 
and a Fe- vacancy ordered BS-AFM insulating phase, re- 
spectively. With this interpretation, the metallic phase 
in the high pressure regime is not necessarily due to the 
charge transfer between two iron sites of different occu- 



pancies as previously expected, but due to the Neel-FM 
phase which respects the ^/5 x ^/5 vacancy ordering. Of 
course, the accompanied structural distortion, in partic- 
ular the difference in zpe, may complicate the comparison 
with experiments. 
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SUPPLEMENTARY INFORMATION 



Definition of magnetic patterns 



We consider several possible ground states with magnetic configurations or patterns as illustrated in FIG. (|S-1[) . and 
we list their total energy under specific pressure in TAB. (|S-I[) . These states can be captured by the extended J1-J2 
Heisenberg model in the classical limit 



<i,j> intra <i,j> inter 



j2SiSj + j2SiSj 



<<i,j>> intra 



<<i,j>>inter 



(a)AFMO 



(b)AFMl 



(c)AFM2 



(d)AFM3 (e)AFM4 (f)AFM5 (g)AFM6 

FIG. S-1: The magnetic structures considered. 



(h)AFM7 



TABLE S-I: Total energies (in eV) per unit cell of different magnetic configurations at different pressures. 



P(GPa) 


AFMO 


AFMl 


AFM2 


AFM4 


AFM5 


AFM6 


AFM7 





-120.163 


-121.152 


-120.709 


-120.517 


-120.559 


-120.719 


-120.583 


2 


-113.772 


-114.625 


-114.201 


-114.076 


-114.094 


-114.175 


-114.094 


4 


-107.688 


-108.397 


-108.010 


-107.956 


-107.932 


-107.942 


-107.919 


6 


-101.841 


-102.405 


-102.061 


-102.088 


-102.013 


-101.956 


-101.994 


8 


-96.207 


-96.607 


-96.313 


-96.431 


-96.295 


-96.178 


-96.289 


10 


-90.767 


-90.980 


-90.750 


-90.955 


-90.766 


-90.581 


-90.798 


12 


-85.489 


-85.510 


-85.494 


-85.635 


-85.496 


-85.138 


-85.505 


14 


-80.366 


-80.169 


-80.364 


-80.450 


-80.364 


-79.834 


-80.366 


16 


-75.355 


-74.955 


-75.355 


-75.526 


-75.355 


-74.685 


-75.355 
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Energy expressions 



The eigen energies of each states per block are listed as follows: 



-Eafmo 
-Eafmi 

-EafM2 

-Eafms 

-E'afM4 

-Eafms 

-^AFMe 

-Eafm? 



(-4Ji - 2J[ + 2J2 + U^)S^ 


(4Ji - 2 J{ + 2 J2 - 4J^)52 
-2J252 

(-4Ji + 2J{ +2J2 -47^)52 

-2J[S^ 

{2J[ - 2J2)S^ 

{-2J[ - 2J2)S^ 



(SI) 
(S2) 
(S3) 
(84) 
(85) 
(S6) 
(87) 
(88) 



Bi-collinear phase cannot be ground state 



In our discussions, the conventional FM phase as weh as the bi-cohinear AFM phase (AFMl, Fig. IS-lf b)) observed 
in some 11-type iron chalcogenides are excluded from the ground state candidates. 

Firstly, we have performed simulations at 12 GPa, confirming that the bi-collinear and FM phases are about ~ 13 
meV/Fe and ~ 19 meV/Fe higher than the Neel-FM phase, respectively. 

Secondly, based on the extended J1-J2 model, it is straightforward to show that the bi-collinear phase could not be 
the ground state when the y/b x ^/5 vacancy superstructure is preserved. 

Proof: If the bi-collinear phase (AFMl) were the ground state, we could obtain from Eqn. [Sl] and Eqn. [S3l 
-4Ji - 2J{ + 2 J2 + 4J2 > and 4Ji - 2J[ + 2 J2 - > 0, respectively. Adding them together yields - J{ + J2 > 0. 
However, this condition contradicts Eqn. [S7l which yields J{ — J2 > 0. Therefore, the bi-collinear phase could not be 
the ground state of the extended J1-J2 model for any Js. The extra next-nearest-neighbor exchanges J3 and J3 are 
needed for this state to be the ground state as in the conventional Ji- J2-J3 model. 



We have performed simulations using 5x5 supercell (230 atoms in all) with 3 randomly distributed Fe-vacancy 
patterns to verify that the vacancy superstructure remains robust under high pressure at least up to 12 GPa. At 
12GPa, the disordered vacancy superstructures are at least 4 meV/Fe higher than non-magnetic >/5 x ^/b phase, and 
are at least 23 meV/Fe higher than the Neel-FM phase. 



For large positive J2 or negative J2, the ground state could be either the BS-AFM state or the Neel-FM state. 
The competition is given by £^bs-afm — £^Neei-FM = £^afm2 — Eafma = 4S'^(2Ji — J[). Hence the critical point is 



In the Neel-FM phase, the band structure indicates an electron pocket for majority spin and a hole pocket for 
minority spin around F. Two large electron pockets can also be identified for the majority spin around X, where no 
structure exists for the minority spin. If the occupancy of Tl sites is reduced by 20% (Tlo.8Fe1.6Se2), the electron 
pocket around F disappears due to the hole doping, and two hole pockets form around F. In this case, one another 
hole pocket also shows up around X. The electron states around Ef are dominated by the Fe-3d orbitals, while the 
Fe-3d and Se-4p orbitals hybridize over a wide range from E^-S eV to Ef^2 eV. 



Stability of vacancy superstructure at high pressure 



Critical point for the BS-AFM/Neel-FM phase transition 



Ji = J[/2. 



Band structure and DOS of the Neel-FM phase 
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FIG. S-2: a) Band structure and b) DOS of Neel-FM TlFei.6Se2 at 12GPa. Notice that majority (solid red line in band 
structure; upper panel in DOS) and minority (dashed blue line in band structure; lower panel in DOS) spins are not degenerate 
in the kx-ky plane, but are degenerate along kz- 



